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STOCHASTIC HODEIXING OF DISPERSION FROM SINGLE ELEVATED SOURCES. £. Robertson 
and P.J. Barry, Atonic Energy of Canada Limited. Chalk River Nuclear 
Laboratories, Chalk River, Ontario, KOJ IJO 

Shore range (e.g. distances <50 ka) atnospherlc dispersion models have 
tradlclonally been used for 'wo sain purposes: a) to predict real time 
concentrations during accidental releases of noxious gases, and b) to predict 
long-term average concentrations for the purpose of controlling pollution 
levels from single sources during routine operations. 

For several decades Gaussian Plume dispersion models together with published 
default dispersion parsmeters have been used for both purposes. However, 
during the last decade an increasing number of studies have shown that the 
reliability of such models and parameter values la poor when used to estimate 
pollutant concentrations for averaging periods of 1 hour or less. 

These studies have been restricted to comparison of predicted to observed 
short-term concentrations in real clme. but, by implication it is often 
Implicitly assumed the model will be unreliable on the longer term as well. 
We have been conducting studies around an isolated stack that continuously 
released a short-lived radioactive noble gaa, Argon-41, at a steady rate. 
This work has confirmed that on a ori-to-one basis in real time the Gaussian 
plume model is a poor predictor for shorter averaging tlmea with discrepancies 
by a factor of 10 or mora occurring 80% of the time. However as the averaging 
time Increases the reliability increases also and when the averaging tine 
exceeded a month the predictions were acceptable. Annual predicted averages 
agreed with those observed within a factor of 2. It was also found that the 
cumulative frequency distributions of predicted hourly averaged concentrations 
agreed reasonably well with tha observed diatrlbutlons. Thus while the model 
gave relatively poor predictions of individual hourly average concentrations 
in real time it nevertheleaa yielded a reasonebly correct ensemble of 
concentrations. 

Thus it follows that, for prediction of long-term average concentrations, or 
the frequency diacrlbuClon of concentrations, the model need not necessarily 
be driven by real time meteorological data, but rather by data synthesized 
from their frequency distributions, %Alch are available In cllmatologlcal 
summaries. Tb« latter method Is easier and thus less costly. 

A stochastic dispersion model based on these Ideas was thersfore developed 
using the following methodology. The code performs Monte-Carlo slmuletions 
with varisble Input meteorological data fed to a standard Gaussian plume 
dispersion equation. The needed meteorological variables are wind speed and 
direction and a measure of the atmospheric stability. The cllmatologlcal 
summary should provide the frequencies of wind directions on es fine an 
angular degree Interval aa possible; however, in the cess of data derived from 
routine synoptic statlona sectors of 22.5" intervals will have to be accepted. 
The cumulative frequenclea are calculated going from 0" to 359" in whatever 
sector Intervals are available. Tha resulting curve is fitted to a polynomial 
equation of the type 
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« + bF + cF" + dF* 



Hh*r« D la th« vlnd dlractlon, F tha cuaulatlv* fraquancy darlvad fron the 
data, and a, b, c -- ate. conatanta. Hoat coaputlng ayatana aupport routines 
for tha evaluation of the conatanta In a polynoalal of any reaaonable degree. 
The polynoalal aquation la written Into tha coda. To aelact a wind direction 
for each -trial-, a randoa nunber froB to 1 la generated, aubatltuted for F 
into the above equation and D calculated. The coaplexlty of the ahape of the 
wind direction frequency curve detcralnea the order of the polynonlal that 
will give a reaaonable approxlutlon to the one obaerved. For the Chalk Blver 
alte, where the wind la atrongly chantwllad Into two relatively narrow angular 
banda, a 9th order waa required. 

A alallar procedure la followed to obtain valuaa for wind apeed and atablllty. 
Stability can be obtained ualng aynoptlc atatlon data and the -STAR- prograa 
froB which the occurrence frequendea of claaaaa A to F can be obtained. 
Stability obtained thle way la a dleerete variable ao that particular class la 
selected froB the randoB nuaber within ranges. For exaBple auppoae classes A, 
B and C occur with frequendea 2%, 10% and 15% reapeetlvely. Then Class A Is' 
selected whenever the randoB nuaber Ilea between and 0.02, claaa B when It 
Ilea between 0.02 and 0.12 and so on. Since wind apeed and atablllty, as 
defined by the STAR prograa, are directly related, t wind apeed distribution 
frequendea are used. 

Once the claaa has been aelectad, valuaa for the horlcontal and vertical 
dlsperalon paraaetera can be calculated for the receptor dlatance froa the 
equatlona provided by Brlgga (1974) or Hoaker (197*). Theae values together 
with a atack aalaalon rata for the hours are iiaed to solve the blvarlate 
Gauaalan pluae equation and obtain a concentration for each receptor. The 
proceas la repeated aa aany tlaea aa neceaaary to obtain a aaooth frequency of 
occurreneea of concentratlona or a arable eatlaate of the long tera aean. 

If the atack aalaalon rate la Independent of waather conditions It too aay be 
selected randoaly froa a frequency dlatrlbutlon of valuaa. 

This aodel la to be applied to other sites to teat Its uaefulness generally. 
Aa a flrat atep It la to be applied to the aaalter stack at Sudbury and the 
predlctlona derived froa It coapered with concentratlona of SO, aeaaured at 
aeveral locatlona in the aurroundlng dlatrlet. 

If the aodel la found to be widely velld. It provldea ^ reletlvely alaple and 
straightforward way of calculating tha stochastic characteristics for a 
cllaatology of elr pollution events. 

Predicting the long tera aean la a uaaful product but knowing tha frequencies 
of occurreneea of different pollution levala about the aean enebles regulators 
to also control tha frequency of occurrence of toxic or nuisance levels. 
Alternatively If air pollution levala are regulated by applying aeteorologlcal 
control over sources, both regulatora and regulated nerd to know the frequency 
and hence the potential coats of production cutbacks aaaoelatad with any given 
strategy for pollution abateaent. 



Brlgga, G.A., 1974, 'Diffusion eatlaatea for aaall ealaalons*, Envlronaental 
Research Laboratorlaa Air Rasourcea, Ataosphere, Turbulence and Diffusion 
Laboratory, Annual Report 1973 USAFC Report ATDL-106, National Oceanic and 
ACaospherlc Adalnlatratlon. 

Uosker, R.P. Jr., 1974, -Eatlaatea of dry deposition and pluae depletion over 
foresta and grasslands', in Froc. of IAEA Syaposlua, -Fhyalcal behaviour of 
radioactive contaalnanta In the etaosphere*, Vlenne, Hoveaber 1973, IAEA, 
Vienna. 



304 



305 



